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Abstract 

We introduce shape fluctuations in a liquid-crystalline system by considering an elementary 
Maier-Saupe lattice model for a mixture of uniaxial and biaxial molecules. Shape variables are 
treated in the annealed (thermalized) limit. We analyze the thermodynamic properties of this 
system in terms of temperature T, concentration c of intrinsically biaxial molecules, and a parameter 
A associated with the degree of biaxiality of the molecules. At the mean-field level, we use standard 
techniques of statistical mechanics to draw global phase diagrams, which are shown to display a 
rich structure, including uniaxial and biaxial nematic phases, a reentrant ordered region, and many 
distinct multicritical points. Also, we use the formalism to write an expansion of the free energy in 
order to make contact with the Landau-de Gennes theory of nematic phase transitions. 
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I. INTRODUCTION 


The characterization of biaxial nematic phases in a number of thermotropic liquid- 
crystalline systems stimulated a revival of interest in the investigation of theoretical 
models to describe biaxial structures [3]. About forty years ago, Freiser |5] showed the 
existence of uniaxial and biaxial nematic phases in a generalization of the mean-field Maier- 
Saupe theory of the nematic transition with the addition of suitably asymmetric degrees of 
freedom. A nematic biaxial phase has also been shown to exist in a lattice model with steric 
interactions between platelets [6], and in a number of calculations for model systems with 
soft and hard-core interactions The early experimental results, however, referred to a 

lyotropic liquid-crystalline mixture PI. which should be better represented by a model of 
uniaxial nematogenic elements ITT, 12], and which motivated the use of an elementary ver¬ 
sion of the Maier-Saupe theory mm to investigate a lattice statistical model for a binary 
mixture of cylinders and disks. We now propose an extension of this elementary model, 
along the lines of Freiser’s generalization of the Maier-Saupe theory, in order to analyze the 
global phase diagram of a mixture of uniaxial and biaxial molecules. 

In some analytical [16] and numerical nn calculations, it has been pointed out that 
shape fluctuations play an important role in the stability of the biaxial nematic phases. In 
spite of the complexity of the liquid-crystalline systems, whose complete description may 
require the introduction of more realistic, and necessarily involved, theoretical models, we 
believe that there is still room for the investigation of elementary statistical lattice models, 
with the addition of some ingredients that may be essential to describe the main features 
of the thermodynamic behavior. Along the lines of Freiser’s early work, we then add extra 
degrees of freedom, of biaxial nature, to an elementary lattice model, which leads to the 
definition of a six-state Maier-Saupe (MS6) model. This MS6 model is similar to an earlier 
proposal by Bocarra and collaborators [7], and may be regarded as a generalization of a 
previously used three-state Potts model to describe the uniaxial nematic transition im. 
Shape fluctuations are taken into account by introducing a “biaxiality parameter” A, and 
by considering a binary mixture of molecules with A = 0 (intrinsically uniaxial molecules) 
and A ^ 0 (intrinsically biaxial molecules). This model system is sufficiently simple to be 
amenable to detailed statistical mechanics calculations for either quenched |l8j or annealed 
mixtures of molecules. We then carry out calculations to obtain global phase diagrams, 
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and write an expansion of the free energy to make contact with the standard form of the 
phenomenological Landau-de Gennes theory of phase transitions. 

This paper is divided as follows. In Section [TTJ we define the MS6 model for a binary 


mixture of molecules, and formulate the statistical problem. In Section III we analyze the 
mean-held equations, draw a number of characteristic phase diagrams, and make contact 


with the Landau-de Gennes theory. Section IV is devoted to a final discussion and to some 
conclusions. 


II. THE SIX-STATE MAIER SAUPE MODEL 


The standard formulations of the Maier-Saupe theory of nematic phase transitions [19] 
can be described in terms of the Hamiltonian 


n = -eJ2 Y, W*’ (D 

(ij) a,0—1,2,3 


where e is a positive parameter, (i,j) means that the sum is over pairs of molecules at sites 
i and j, and is the symmetric traceless quadrupole tensor associated with a molecule at 
site i. In general, fl, may be written in terms of the direction cosines or Euler angles which 
connect the laboratory and molecular frames. From the traceless condition, we write the 
eigenvalues Ax = —1 +A, A 2 = —1 — A, and A 3 = 2, of the tensor G, where A is a parameter 
that gauges the degree of biaxiality. 

The problem is considerably simplified if we resort to a discretization of directions, which 
has been used to describe the isotropic-nematic transition pn. We then assume that the 
principal molecular axes are restricted to the directions of the Cartesian coordinates of the 
laboratory. Therefore, the quadrupole tensor Q can assume only six states represented by 
the matrices 
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which leads to the definition of the six-state Maier-Saupe (MS6) model. If the molecules are 
intrinsically uniaxial (A = 0), we regain a three-state model, which has been used to describe 
the transition from the isotropic to the uniaxial nematic phase ra. and to investigate the 
existence of a biaxial nematic phase in a binary mixture of cylinders and disks [T5HT5] . 

The thermodynamic behavior of the MS6 model is determined from the canonical parti¬ 
tion function 


Z = E exP 

(fid 


^E E 

(ij) « =1 -2,3 


aa 
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(3) 


where /3e = 1/T, so that T is the temperature in suitable units, and the first sum is over all 
microscopic configurations {fij} of this six-state model. This problem is further simplified 
if we consider a fully connected model, with equal interactions between all pairs of sites. At 
this mean-field level, if A = 0, we anticipate just a first-order transition between an isotropic 
and a uniaxial nematic phase. If A ^ 0, however, we can describe the transition to a stable 
biaxial nematic phase. 

We now turn to a mixture of intrinsically uniaxial (A = 0) and intrinsically biaxial 
(A 7 ^ 0) molecules. In this mixture we have two sets of degrees of freedom: (i) orienta¬ 
tional degrees of freedom, {fij}, of quadrupolar nature, and (ii) shape-disordered degrees of 
freedom, {Aj}, with either Aj = 0 or Aj = A ^ 0, at all lattice sites. These two sets of 
degrees of freedom may be associated with quite different relaxation times, which leads to 
the distinction between annealed and quenched situations PSIEDJ. In the quenched case, 
the “shape-disordered” degrees of freedom never reach thermal equilibrium during the exper¬ 
imental times. Given a configuration {Aj}, we calculate a partition function Z = Z ({Aj}), 
and a configuration-dependent free energy, / ({Aj}). The free energy of the system is an 
average of / ({Aj}) over the shape-disordered degrees of freedom, and the concentration 
of intrinsically biaxial molecules is not a true variable of equilibrium thermodynamics. In 
the annealed case, the two sets of degrees of freedom are supposed to thermalize during 
the experimental time, so that concentration and chemical potential are thermodynamically 
conjugate variables. In the annealed case, given the concentration, both types of particles 
are free to move across the system in order to minimize the free energy. In this work, we 
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consider annealed disorder only, which is more appropriate to a liquid-crystalline system. 

In the annealed case, consider a binary mixture of N { intrinsically biaxial molecules 
(A 7^ 0) and N 2 = N — Ni uniaxial molecules (A = 0). Given the numbers of uniaxial and 
biaxial molecules, the canonical partition function is a sum over orientational and disorder 
configurations, 


Za = 


{«<}{ A d 


AE (a<) si-( a, 

(ij) a=l,2,3 


(4) 


where the prime in the second sum indicates the restriction 

N 

E^-V.A (5) 

2—1 

At this stage, it is convenient to introduce a chemical potential /j and change to a grand 
ensemble. First, we redefine the shape variable of molecule % such that 


A i = mA, 


( 6 ) 


where 

0, for a uniaxial object, 

1, for a biaxial object. 
Then, the grand partition function is given by 


Hi = 


(7) 


EE 
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where /i is the chemical potential that controls the number of biaxial molecules. We now 
remark that the sums over configurations in Eq. (|8]) are no longer restricted, which makes 
it possible to carry out the calculations in the mean-held limit, as it is detailed in the next 
Section. 


III. MEAN-FIELD CALCULATIONS 


The mean-held version of the MS6 model is given by the Hamiltonian 
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The grand partition function 
identities, 


exp 


in Eq. ([ 8 | can be factorized by using three Gaussian 
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with a G {1,2,3}. This factorization effectively decouples the problem of calculating the 
grand partition function, and the sums over and {rii} can be performed in a straight¬ 
forward way, so that we can write 
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where ifi a functional of {q a }. In the thermodynamic limit, the integral can be calculated 
by standard saddle-point techniques. Thermodynamic equilibrium is then associated with 
the minimization of ifi with respect to {q a }, from which we obtain self-consistent mean-field 
equations for these quantities. These equations show that q\ + q 2 + <?3 = 0, which suggests 
the introduction of a symmetric traceless tensor, 


^ gi 0 0 ^ 
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(13) 


0 q 2 0 
y 0 0 q 3 

as an appropriate thermal average of fij(nj). Using the traceless condition, it is convenient 
to rewrite Q as 

/ -S-7] 0 0 

0 -S + 7j 0 
0 0 2S 

in terms of two scalar parameters, S and rj. The isotropic phase is given by S — rj — 0. The 
nematic uniaxial phase is given by S' 7 ^ 0 and r] = 0 (or ij = ±3S). In the biaxial phase, we 
have S 7 ^ 0 and 7] 7 ^ 0. 

In the following paragraphs we write explicit expressions for the thermodynamic poten¬ 
tials of the uniform system and of the annealed binary mixture. From these expressions, 
it is easy to perform numerical calculations to draw a plethora of global phase diagrams in 
terms of the model parameters. In order to asymptotically check the numerical Endings, and 
to make contact with established phenomenological results, we may also write an expansion 
of the thermodynamic potential in terms of the invariants of the tensor order parameter, 
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I n = TrQ n , with n — 1, 2, .... Due to the symmetry properties of Q, all these invariants 
can be written as polynomials depending on two basic invariants, given by 


/ 2 = Tr Q 2 = i (3S 2 + if), 


(14) 


and 

I 3 = TvQ 3 =^S(S 2 -7f). 

Therefore, the usual form of the Landau-de Gennes expansion is written as 


A T B T C 2 D E 3 E' 2 

9 - 9a + ^ h + j h + j h + -jr hh + jr h + -jr 4 + - 


(15) 


(16) 


According to this phenomenological expansion [2T|, there is a Landau multicritical point for 
A = B = 0. In the vicinity of this Landau point, we can establish parametric expressions 
for the lines of phase transitions between the isotropic and nematic phases. 

We now consider the specific cases of uniform and annealed systems. 


A. Uniform case 


Consider a system of intrinsically biaxial molecules (Aj = A for all i). Assuming the 
discretization of orientations, and setting /i = 0 in Eq. (J 8 ]) , since disorder plays no role, the 
functional -0 is written as 
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where T = (0e) -1 . Minimizing 0 with respect to S and 77 leads to the self-consistent mean- 
field equations, S = F\ (S, ?/; T, A) and r] = F 2 (S, 77 ; T, A). The values of S and r) at the 
absolute minimum of 0 correspond to the thermodynamic equilibrium values for a fixed 
temperature T and degree of biaxiality A. The free energy f = f (T, A) of the system is 
obtained from 0 by inserting the equilibrium values of S and 77 . 

Figure [l] shows the phase diagram in the T — A plane, which is obtained by solving 
the mean-field equations numerically. As it should be anticipated from phenomenological 
arguments, this phase diagram shows two lines of continuous transitions (solid lines) from a 
biaxial nematic region to the N + (prolate) and N_ (oblate) uniaxial nematic regions. These 
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FIG. 1: Phase diagram, in terms of temperature T and the degree of biaxiality A, for a system 
of intrinsically biaxial molecules. N + and N_ are uniaxial nematic prolate and oblate phases, 
respectively. B is the nematic biaxial phase, L is the Landau multicritical point, and ISO is the 
isotropic region. Continuous lines indicate continuous transitions and dashed lines represent first- 
order phase transitions. 

critical lines meet at a Landau multicritical point, L, on the first-order boundary (dashed 
lines) between the isotropic and the uniaxial nematic phases. It should be remarked that we 
regain an intrinsically uniaxial system for A = 3. The phase diagram for A > 3 is mapped 
onto the region 1 < A' < 3 by the transformations A' = (A + 3) / (A — 1) and T 1 = 
4 Tj (A — l) 2 . Note that a similar model for asymmetric ellipsoids leads to essentially the 
same type of phase diagram [Tj . Also, a number of calculations for continuous orientational 
degrees of freedom lead to the same characteristic topology of this phase diagram (see, for 
example, the works of Luckhurst and collaborators |9] and of Xheng and Palffy-Murhoray 


From the expression of the free energy, we obtain the parameter-dependent coefficients 
of a Landau-de Gennes expansion about the Landau multicritical point, 


A = 1 


3 +A 2 _ 
T ’ 


S = 9/A2 
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J 1 '! 


T 3 \ 4 


C = — \ - + -A + -A 4 
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( 20 ) 


D = ——- (A 4 + 2A 2 - 3) , 

16T 4 v ’ ’ 

E = - 48( ' )y . 5 (41A 6 + 315A 4 + 1215A 2 + 1053) , (21) 

E' = —(A 6 + 225A 4 - 405A 2 + 243) . (22) 

Therefore, the Landau point L is located at Tj, = 4 and Al = 1. In the vicinity of L, limiting 
to first order terms in (T — 4) and (A — 1), we have A = (1/4) (T — 4) — (1/2) (A — 1), 
B = (9/16) (A - 1), C = 1/16, D = 0, E = -41/7680, and E' = 1/640. The sign of 
E' indicates the stability of the biaxial nematic phase near the Landau point [19, I2T] , In 
this mean-field scenario, at fixed A / 1, as the temperature decreases from a sufficiently 
large value, the system goes from an isotropic phase to a uniaxial nematic phase, and then 
to a biaxial nematic phase, according to the prediction of the early work of Freiser [5]. It 
should be remarked that the phenomenological Landau parameters are written in terms of 
the parameters of the underlying molecular model, which makes it easier to investigate a 
large range of values. 


B. Annealed disorder 

In the annealed case, we calculate 0, given by Eq. in terms of the proper ther¬ 

modynamic field variables, temperature T and chemical potential // = eTlnz, where z is 
a fugacity. For the fully-connected MS 6 model of a binary mixture of uniaxial and biaxial 
molecules, the functional ijj is given by 

-0 = - (3 S 2 + r/ 2 ) — T In 2 — T In cr, (23) 

where 


_3 (S+rj) 

a = e °- T <1 + z cosh 
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(24) 


Again, the thermodynamic stable values of S and r] are chosen to minimize the function 0 
for fixed values of temperature T, degree of biaxiality A, and chemical potential //. 
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Inserting the equilibrium values of S and 77 into the expression of '0, we obtain the grand 
potential as a function of T, /i, and A. If we wish to work with a fixed concentration c of 
the intrinsically biaxial molecules, the free energy 0 (c, T; A) comes from the definition 

0 (c, T; A) = 0 + clnz, (25) 


where the fugacity is eliminated by the expression 


c = 


<90 


(26) 


and we should insert the equilibrium values of S and rj (with the proviso of a Maxwell 
construction whenever it is necessary). 

The grand potential can be used to write a Landau-de Gennes expansion with coefficients 


I , 3 + ^(3 + A 2 ) 9 2 (A 2 — 1) — 1 

T(l + z) ’ 2 T 2 (l + z) ’ 
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where 

d 0 = 486 - 405A 2 - 45A J + A 8 . 
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(34) 
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and 

d x = 243 - 405A 2 + 225A 4 + A 6 . (36) 


From the usual condition A = B = 0, we locate a Landau multicritical point, 


dL = ~Tl In (A 2 — l) , T l = 4, 


(37) 
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FIG. 2: Energy levels as a function of the degree of biaxiality A. Level crossings account for 
variations in the entropic contribution to the free energy. 



A 

which requires A 2 > 1, and gives an indication of the existence of qualitatively different 
phase diagrams (as there is no Landau point for A 2 <1). 

It is easy to show that this annealed version of the binary mixture exhibits phase diagrams 
with many distinct topologies. This can be anticipated from an analysis of the energy levels 
associated with the interaction between molecules, as shown in Figure [2j In fact, these 
energy levels are associated with different degrees of degeneracy oj, and some levels cross 
each other as the biaxiality of the molecules is changed. These degeneracies account for 
entropic contributions to the free energy, which do affect the equilibrium phase behavior 
of the system. Therefore, we anticipate qualitative changes in the phase diagrams as A 
assumes values close to the location of the energy level crossings. 

We now discuss the various topologies exhibited by the phase diagrams as the biaxiality 
parameter is changed. 

Figures [3]|4] show T — n and T — c phase diagrams for fixed degree of biaxiality, A = 1, 
where /j is the chemical potential and c is the concentration of intrinsically biaxial molecules. 
From the thermodynamic point of view, the first-order boundaries (dashed lines) in the T—fi 
plane are mapped into coexistence regions (gray regions) in the T — c plane. For low temper¬ 
atures and intermediate concentrations, there is a coexistence region between the uniaxial 
prolate (N + ) and biaxial (B) phases. However, at higher concentrations and intermediate 
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FIG. 3: Phase diagram for fixed degree of biaxiality A = 1, where T is temperature and fi is 
chemical potential. There is a tricritical point, TC, along the boundary separating biaxial (B) and 
uniaxial prolate nematic (N_|_) phases. There is no direct phase transition between the isotropic 
(ISO) and biaxial phases (for finite values of the chemical potential). 

temperatures, there is a second-order phase transition (continuous line) between the N + and 
B phases. In fact, the phase diagram exhibits a tricritical point, TC, along the boundary 
between N + and B. At higher temperatures, there is a first-order phase transition between 
N + and the isotropic (ISO) phases, with a very thin coexistence region, as shown in the in¬ 
set. Note that an incipient Landau point appears at c = 1, which corresponds to an infinite 
chemical potential, in agreement with the Landau-de Gennes expansion. Phase diagrams 
with a similar topology (but with no Landau point) can be drawn for 0 < A < 1. 

Figure [5] shows the phase diagram for a fixed degree of biaxiality A = 1.4. Similar to 
Figure[4| the N + and B phases coexist for intermediate concentrations and low temperatures. 
However, the system also displays a uniaxial oblate (N_) nematic phase, which appears at 
higher concentrations and intermediate temperatures. The B phase appears between the two 
uniaxial phases, at intermediate temperatures. All three ordered phases become identical 
to the ISO phase at a Landau multicritical point, L. Also, note that the B phase presents a 
discrete reentrant behavior close to L. Note that the changes in the topology of the phase 
diagrams shown in Figures [4] and [5] are in agreement with the dependence of the energy 
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FIG. 4: Phase diagram in terms of the temperature (T) versus the concentraton (c) of biaxial 
molecules, at a fixed degree of biaxiality, A = 1. There is an incipient Landau multicritical point at 
c = 1, where the two ordered phases and the isotropic phase become identical. Uniaxial (N+) and 
biaxial (B) phases coexist at low temperatures and intermediate concentrations (gray region). The 
system exhibits a tricritical point, TC, along the boundary of the B phase. The inset represents a 
zoom of the coexistence (gray) region between the uniaxial nematic and isotropic phases. 

levels on the degree of biaxiality A (see Figure [2]). 

The phase behavior of the system changes significantly for degree of biaxiality around 
A = 1.5, which is close to another crossing of energy levels. For example, in Figure[6]we show 
the phase diagram for degree of biaxiality A = 1.54. The low temperature biaxial phase B_ 
is represented by a tensor order parameter Q whose largest eigenvalue (in absolute value) 
is negative. However, an additional biaxial phase B appears near the Landau point. The 
two biaxial phases are stable in disconnected regions of the phase diagram. There is then a 
coexistence region between the uniaxial nematic phases N + and N_. This region is limited 
by two critical end points, CE, associated with the biaxial phases. Also, the system exhibits 
a tricritical point TC related to the coexistence region of N + and B structures. As in the 
case A = 1.4, shown in Figure [5j the biaxial phase is reentrant in the vicinity of the Landau 
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FIG. 5: Phase diagram, in terms of temperature T and concentration c of biaxial molecules, for 
degree of biaxiality A = 1.4. There is a Landau point L. Note that the biaxial phase displays a 
reentrant behavior near L. 

point. 

The topological changes in the phase diagram can be represented by the projections of 
the lines of different mnlticritical points on the T — A plane, as indicated in Figure [8j The 
temperature of the Landau point (black) is a constant function of the degree of biaxiality 
A, which is in agreement with a Landau expansion. The temperature of the tricritical point 
(green) increases monotonically with A, as suggested by Figures [4j|6] However, the line of 
critical end points (blue) presents a reentrant behavior in the vicinity of A = 1.55, giving 
rise to the two critical end points, as it is shown in Figure |6j All multicritical lines meet 
at a higher-order multicritical point Mi. Note that the tricritical and high-temperature 
critical end points are associated with the biaxial phase near the Landau point. As a result, 
there is no stable biaxial phase in the vicinity of L for A > A m 1 — 1.5576, although the 
low-temperature biaxial phase survives in a small region at high concentrations and low 
temperatures. This is illustrated in Figure [7J for A = 1.57. There is no stable biaxial 
nematic phase close to the Landau point L, which marks the meeting of various coexistence 
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FIG. 6: Phase diagram in the c — T plane for degree of biaxiality A = 1.54, where T is the 
temperature and c is the concentation of biaxial objects. Two nematic biaxial phases are stable 
in disconnected regions of the phase diagram. There is a coexistence region between the nematic 
uniaxial phases, ending at two critical end points CE. The dotted horizontal lines represent special 
tie lines associated with the critical end points. 

lines separating the isotropic phase and the two nematic uniaxial phases. Along those lines, 
there is a coexistence of phases with different values of nematic order parameter S, but the 
size of the coexistence region tends to vanish as we approach the Landau point. 

As the degree of biaxiality is further increased, the system exhibits other distinct phase 
diagrams, as seen in Figures |9]fT0| According to Figure |2j multiple level crossings occur at 
A = 3, which is a strong suggestion of changes in the phase behavior of the system. 

Figure [9] shows the phase diagram for A = 2.9. There appears a triple point associated 
with the coexistence of two uniaxial nematic oblate phases and a uniaxial prolate phase. 
Also, the coexistence region of uniaxial nematic oblate phases ends at a simple critical point 
C. Although it is not shown in Figure [9j a stable biaxial nematic phase is still present, at 
low temperatures and high concentrations, as well as the critical end point associated with 
the biaxial phase. 

As A is further increased, the simple critical point moves upward in the phase diagram, 
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FIG. 7: Phase diagram in terms of temperature (T) versus concentrarion (c) of biaxial molecules 
for fixed degree of biaxiality, A = 1.57. The biaxial phase B_ is located at low temperatures and 
high concentations. There is no stable B phase near the Landau point L. The nematic uniaxial 
phases coexist in a region ending at a critical end point CE. The dotted horizontal line represents 
a special tie line associated with the critical end point. 


approaching the lower border of the coexistence region between the uniaxial N_ and ISO 
phases. Then, the simple critical point is replaced by a second triple point. For example, 


Figure 10 represents a phase diagram for A = 3. Due to the second triple point, there 
appears an isolated second uniaxial oblate phase Note that, according to Eq. the 
value A = 3 corresponds to a binary mixture of rods and plates with asymmetric interaction 
energies. A symmetric choice of interaction energies na presents a much simpler phase 
diagram, with single nematic uniaxial prolate and oblate phases, in addition to an isotropic 
phase, and no triple points. 

Figures 11 and [l2| represent the phase diagrams for A = 5 and A = 6, respectively. In 
these diagrams the values of concentrations are conveniently rescaled, so that visual effects 
are improved. Also, note that we introduce some separations just to emphasize the more 
interesting sectors of this phase diagram. There is no stable Landau point in both phase 


diagrams. In Figure 11, the uniaxial phases N + and N_ coexist with the ISO phase at a 
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FIG. 8: Lines of multicritical points in the A — T plane. There is a higher-order multicritical point 
Mi, at which the lines of Landau points (black), tricritical points (green) and critical end points 
(blue) meet. There is then no stable biaxial phase for A > Am, in the vicinity of the Landau point. 


triple point, and a biaxial B phase remains stable at high concentrations. However, in 


Figure 12, there is a coexistence region of B_ and ISO phases, as well as a triple point of 


coexistence of ISO, B and N + phases. Figures 11 and 12 suggest that the triple points 
meet the critical end point as the degree of biaxiality is increased. 

In Figure [13] we draw the projections of distinct multicritical points on the A — T plane. 
The line of Landau points is still present, as well as the low-temperature part of the line 
of critical end points associated with B . Furthermore, there are lines of triple points 
associated with the ISO and various nematic uniaxial phases, as depicted in Figures [9] and 


10 The line of Landau points meets the triple lines at another special multicritical point, 


which we call a Landau end point, LE. This point is characterized by the coexistence of a 
critical ISO phase and a non-critical N_ phase. Otherwise, two triple lines meet the line of 
simple critical points at two multicritical end points, MCE. In addition, for A ~ 5.5, the 
triple line meets the critical end line at a multicritical point, M 2 , where the critical phase 
N_ coexists with the non-critical phases N + and ISO. Consequently, for A > 5.5, there is a 
coexistence region between ISO and B phases as we increase the concentration of biaxial 
molecules. 
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FIG. 9: Phase diagram in terms of temperature T and concentration c of biaxial objects, for degree 
of biaxiality A = 2.9. There is no biaxial nematic phase around the Landau point L. C indicates a 
simple critical point associated with a coexistence region between two uniaxial oblate phases. The 
dotted line indicates a triple point, corresponding to the coexistence of two uniaxial oblate and one 
uniaxial prolate phases. 

IV. CONCLUSIONS 

We introduced an elementary six-state Maier-Saupe (MS6) lattice model, which is ob¬ 
tained by the addition of extra degrees of freedom, of biaxial nature, to an earlier three-state 
model. We then described a phase diagram with biaxial as well uniaxial nematic structures 
and an isotropic phase. A fully connected MS6 lattice model, of mean-held character, is 
sufficiently simple to be amenable to a detailed treatment by standard statistical mechan¬ 
ics techniques. Results are obtained in terms of a parameter A of shape fluctuations that 
gauges the degree of biaxiality. We then used this MS6 model to consider a binary mixture 
of intrinsically uniaxial (A = 0) and intrinsically biaxial (A ^ 0) molecules, and investigated 
the effects of “shape fluctuations" on the phase diagrams in terms of temperature and either 
chemical potential or concentration of biaxial molecules. 

Taking into account the fluid character of the liquid-crystalline systems, in the present 
work we restricted the analysis to a thermalized (annealed) situation, in which case ori- 
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FIG. 10: Phase diagram in terms of temperature T and concentration c of biaxial objects, for degree 
of biaxiality A = 3.0. This system corresponds to a mixture of rods and plates with asymmetric 
interaction energies. A uniaxial nematic prolate phase, N+, and two uniaxial nematic oblate phases, 
Nl and NH, are present. There are two triple points indicated by the dotted tie lines. 

entational and shape degrees of freedom reach equilibrium simultaneously. We obtained a 
wealth of topologically distinct phase diagrams, with several nematically ordered structures 
and multicritical points. 

In the uniform case, in terms of temperature T and the parameter A, we regained the 
first-order transitions between isotropic and uniaxial nematic phases, and the critical lines 
between the biaxial and uniaxial nematic phases, which meet at a Landau multicritical 
point. For a binary mixture of biaxial and uniaxial molecules, we have drawn a number of 
phase diagrams in terms of temperature versus concentration of biaxial molecules, with fixed 
values of A, which display many distinct features. Depending on parameters, there appear 
additional multicritical points, such as tricritical and critical end points associated with the 
biaxial nematic phase. The Landau, tricritical and critical end points may give rise to a 
higher-order multicritical point. Also, depending on the range of values of A, there may be 
a line of Landau points meeting a line of triple points at a Landau end point. This intricate 
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FIG. 11: Phase diagram in terms of temperature T and concentration c of biaxial molecules, for 
degree of biaxiality A = 5.0. The Landau point is absent. There appears a triple point, at which 
N+, N_ and ISO phases coexist. The B_ phase is stable at high concentations. 
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FIG. 12: Phase diagram in terms of temperature T and concentration c of biaxial objects, for 
degree of biaxiality A = 6.0. The Landau point is absent. There is a coexistence region between 
ISO and B phases. In addition, ISO, B„ and N + phases coexist at a triple point. The B_ phase 
is stable at high concentations. 
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FIG. 13: Lines of phase transitions in the A — T plane. The dashed black lines represent lines of 
triple points. The line of simple critical points (red) ends at two multicritical end points MCE. The 
line of Landau poins meets the lines of triple points at a Landau end point LE. Also, the line of 
critical end points (blue) crosses the line of triples points at a multicritical point M 2 . 

behavior of the mixtures of intrinsically uniaxial and biaxial molecules can be understood 
in terms of crossings of the microscopic energy levels as we change the degree of biaxiality 
A. 

The explicit expressions for the free energy were used to obtain the coefficients of an 
expansion at high temperatures, in the vicinity of the Landau multicritical point, and to 
make contact with the Landau-de Gennes theory. We were then able to check our numer¬ 
ical findings against a number of phenomenological calculations of the literature. Also, we 
provided a simple way of obtaining the expansion coefficients in terms of the values of the 
molecular parameters. 

The calculations predict the presence of stable biaxial nematic phase at low temperatures 
and sufficiently high concentrations of biaxial molecules, which is in agreement with recent 
calculations of Longa and coworkers [T5] for a more elaborate model system of a mixture of 
biaxial molecules in the annealed situation. Depending on the degree of biaxiality, we predict 
first-order transitions between biaxial and uniaxial nematic phases, as well as tricritical 
points, even in the absence of a Landau multicritical point. For larger values of the degree of 
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biaxiality, the Landau point splits into lines of triple points. Also, we note a clear reentrance 
of the biaxial regions for some choices of the parameters, which is agreement with the early 
work of Alben for a lattice model of platelets [6]. 

In some recent publications, Akpinar, Reis, and Figueiredo-Neto |23] |24] reported X-ray 
and optical characterizations of biaxial nematic structures in a large class of quaternary 
liquid-crystalline mixtures. From these measurements, it has been possible to establish 
many novel phase diagrams in terms of temperature and molar fraction of the components, 
which represents a real advance with respect to the early work of Yu and Saupe for a ternary 
lyotropic mixture. For all concentrations of the amphiphile component, if there is a nematic 
biaxial structure, it is thermodynamically stable at intermediate temperatures, between 
regions of different uniaxial nematic structures, at lower and higher temperatures. Also 
there are examples of temperature-concentration phase diagrams with a clear indication of 
the existence of a Landau multicritical point. With a suitable choice of parameters, the MS6 
model of a binary mixtures can qualitatively reproduce all of these observations. In a very 
recent experimental investigation, Amaral and coworkers [22] reanalyzed the phase diagram 
of a ternary SDS lyotropic mixture, and pointed out a peculiar coexistence of uniaxial and 
biaxial nematic structures, which still seems to demand a theoretical explanation. 

The present calculations, for the elementary MSG lattice model, at the mean-held level, are 
a contribution to the understanding of the effects of shape fluctuations on the thermodynamic 
behavior of complex liquid-crystalline systems. The model of a binary mixture is sufficiently 
simple to produce a number of analytical and numerical results for a wide range of values 
of the molecular parameters. Also, it seems to be possible to go beyond the mean-held 
scenario. The use of more powerful techniques may uncover additional aspects of the phase 
diagrams, in particular limitations of the mean-held approach at very low and very high 
concentrations. 
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